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We consider a slowly rotating rectangular billiard with moving boundaries and use 
canonical perturbation theory to describe the dynamics of a billiard particle. In the 
process of slow evolution certain resonance conditions can be satisfied. Correspond- 

(N. 

\ iiigly, phenomena of scattering on a resonance and capture into a resonance happen 

CN 

I in the system. These phenomena lead to destruction of adiabatic invariance and to 
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unlimited acceleration of the particle. 

When one slowly changes parameters of an integrable Hamiltonian system with 
two degrees of freedom, two classical actions of the unperturbed system are ap- 

X! . 

■ proximately conserved as adiabatic invariants. However, if during slow evolution 



resonance relations between frequencies of the system are satisfied, adiabatic invari- 
ance is destroyed. The paper describes in detail how it happens in a rectangular 
billiard system with moving boundaries and how it can lead to unlimited accelera- 
tion (so-called Fermi acceleration) of the billiard particle. 
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INTRODUCTION 



Billiards are important models in the theory of dynamical systems and its applications 



]J]). Recently 
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(Ref. [l-3|), especially in optics and matter- wave optics (see, e.g., Ref. 
billiards with varying parameters became the object of studies ( see, e.g., Ref. 
One of the most interesting questions in this area of research is Fermi acceleration (Ref. 

3)- 

In the present paper we consider acceleration of a billiard particle due to scattering on 
resonances in a rectangular billiard with slowly varying parameters. We use the methods 
of analysis of these resonant phenomena designed for Hamiltonian systems with slow and 



fast variables (Ref. [24|). The methods were developed and the corresponding theorems 



were proved for smooth systems. However, previous studies (Ref. Q, show that 

they can be also applied adequately for billiard systems, which possess discontinuities 
(impacts). In Ref. 16|, acceleration of a particle due to scattering on resonances and 
capture into a resonance was considered mostly for the case of crossing periodically a 
single low-order resonance. The numerical calculations revealed very slow acceleration, 
and a possibility of unlimited acceleration remained an open question. Here we describe 
a geometric obstacle for acceleration that considerably slows it down, so that it is hard 
to detect the energy growth numerically. When this obstacle is not present, the particle 
accelerates without bounds fast enough to clearly see the energy growth in numerical 
calculations. 



II. THE MODEL 



Let us consider a particle in a rectangle 2D box rotating at a constant angular velocity 
a; > 0. We assume that the rotation is slow: a; ^ 1. The Hamiltonian of the system in 
the rotating coordinate frame is 
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FIG. 1: Phase portraits of the Hamiltonian Fq in ([8]). a) \a\ < |6|, the phase portrait has an 
osciUatory domain, b) \a\ > \b\, there is no oscillatory domain in the phase portrait. In both 
plots, a < 0. 
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FIG. 2: Jumps of adiabatic invariant Ii. Parameters of the system: e = 1- 10~'^, uj = 4.5 • 10~^, 
di = dio{l + Ai cos(et)), d2 = (i2o(l + ^2 cos(rf)), where Ai = —0.1, A2 = 0.3, dio = 1, = 1.4. 
Inset: a single jump of the adiabatic invariant. 



H = 2(^1 +P2) +^(Pi92 -P2gi) + U{qi,q2,et), 



(1) 
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where Pi,qi are canonically conjugate momenta and coordinates, respectively, and 
U{qi,q2,et) is the hard wall potential of the 2D box, and £ is a small positive parame- 
ter (the time-dependence arises from motion of boundaries described below). The second 
term in Eq. ([1]) comes from switching to the rotating frame (where billiard is irrotational): 
it is —uL, where L = —{piq2 — P2qi) is the angular momentum of the particle about the 



coordinate origin 



28| . The rotation is needed to couple two degrees of freedom, otherwise 



the system would decouple on two one-dimensional billiards. 



16| and is given here 



The Hamiltonian formalism outlined below is available in Ref. 
for self-consistency of the presentation. Using a canonical transformation we introduce 
new variables (Jj, i = 1, 2 such that 

_ f -di + 2di(j)i/7i, if < 0i < TT, 

~ \3di- 2di(f)i/7r, if vr < 0^ < 27r, ^ ' 

Pi = ^sgn(sin0i), 

where 2di are lengths of the sides of the 2D box. We assume that these lengths slowly 
depend on time: di = di{et), where e ~ w ^ 1. If di are constant, and the box does not 
rotate, the new variables (Jj, (pi) are just the action-angle variables of the system. 
In the new variables the Hamiltonian of the system is: 



8[dl^dl) 



{hdlh + I2dlk2) . , , , . 



{hdlki - I2dlk2) . , . , N 

+ d,d2klkl ""^^^'^^ + ^^'^^^ 



+ e[E^{<Pi) + E2{<P2)] = H^ + ooH^, (3) 



where primes denote summation over odd k. We use Fourier series expansion of expres- 
sions Eq. ([2]) for Pi,qi (see Appendix). Terms Ei appear because the transformation ([2]) 
is non-autonomous; the explicit form of these terms is not important for our discussion. 

Let m, n be relatively prime integers. Let us call the terms in ([3]) with ki = rm, k2 = rn 
and integer r (m, n)-resonant, if the following resonance condition is met: 



5 



Physically, the resonance condition means that frequencies of oscillation of a parti- 
cle between pairs of opposite walls are commensurable {muji — nuj2, where — ^ — 
j — 1,2 are frequencies of oscillation "along" and "across" the rectangular bil- 
liard, correspondingly). If the resonance condition is met in the unperturbed billiard, the 
particle is moving along a closed trajectory. 

The resonance condition defines a corresponding resonance ray in the plane of actions: 

where am-.n is a slow function of time. 



III. PASSAGE THROUGH A RESONANCE 

If a phase point has such values of /i = Ii'^ and I2 — I2'' that in the plane of action 
it is located far from any low-order resonance ray, then in the process of evolution values 
of Ii and I2 will be approximate adiabatic invariants and the phase point will undergo 
small oscillations around (/^'^^ /2°Y 

However, as the time grows, the values of the lengths di change and the system may 
approach a resonance. There, actions are no longer well-conserved, and two phenomena 
can happen: scattering on resonance and capture into the resonance. They are illustrated 
in Fig. 2 and Fig. 3. Fig. 2 presents evolution of action 7i in a billiard with periodically 
oscillating walls. The action undergoes small oscillations remaining almost constant until 
a resonance condition is fulfilled in the system. At a resonance, it may undergo a jump, 
which can be seen in detail in the inset of Fig. 2. This is scattering on a resonance. An- 
other option for a phase point is capture into the resonance shown in Fig.3. As the system 
approaches a resonance, the phase point can be captured in it and continue its motion in 
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FIG. 3: (a) Jumps of adiabatic invariant Ii (due to scattering on a resonance) and a capture 
of the phase point into the 1:1 resonance. The captured point moves in the plane of actions 
with the resonant ray I2 = cti-ili along the resonant trajectory J = const until it escapes 
from the resonance, (b) Behaviour of /i,/2 near and at the capture into the resonance 1:1 
(arcs in the center of the figure correspond to the phase point captured into the resonance) (c) 
Behaviour of the winding number {u;i/uj2) as a function of time. For the point captured into the 
resonance, this winding number stays approximately constant and undergoes small oscillations 
around the resonance value (a;i/a;2=l in this particular case), (d) The same as in panel (c), but 
on a different scale. Scattering on a resonance 1:1 causes a jump in the winding number due 
to jumps in /i,/2 (inset shows a single jump with even greater resolution). Parameters of the 
system: e = 1 ■ 10~^, a; = 4.5 • 10~^, di = dio{l + Ai cos(et)), d2 = (^20(1 + vl2Cos(et)), where 
Ai = -0.1, A2 = 0.28, dio = 1, d2o = 1.4. 
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FIG. 4: Diffusion in the action space (/i, 12) due to periodic crossing of resonance lines (schemat- 
ically), (a) Only one low-order resonance is crossed. The particle can diffuse in the action space 
due to kicks it receives at each passage of the resonance ray IJ*^* = am-n{t)Ii through it, but dif- 
fusion happen only along an interval of the resonance trajectory J =const within the resonance 
sector since all these kicks are parallel to each other and are restricted by the resonance sector. 
Solid lines shows the resonance rays I2 = am-nh for two values of the slope am,n, that is, amm 
and amax- These two lines define the resonance sector. Dashed lines illustrate the resonance 
sector of another resonance, which our particle does not have chance to cross, (b) Resonance 
sectors of several low-order resonances overlap, correspondingly the particle consequently cross- 
es several resonances. At crossing of either resonance the phase point receives a kick along the 
direction corresponding to the resonance. Since these directions are not parallel (J = nli + m/2 
is defined differently at different resonances), unlimited diffusion is now possible. Note that each 
resonance is crossed as a singe resonance, 'overlapping' means that some sector in the action 
plane is sweeped by several resonance rays, but in different times. 

such a way that the resonance condition would persist. This is clearly demonstrated in 
Fig. 3c, where a winding number Ui/u2 is shown. Capture into the 1 : 1 resonance leads 
to a remarkable conservation of the winding number during a long time (comparable to 
a slow period 2tc /e). There are also several events of scattering on a resonance which are 
clearly seen in Fig.3a (as jumps in the actions similar to Fig. 2), and can also be noticed in 
behaviour of the winding number in Figs. 3c, d. Physically, captured point moves in such 
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FIG. 5: Diffusion in the action space, (a) Only one low-order resonance influences dynamics of 
the particle, crossings of the higher-order resonances nearby lead to slight perturbation of the 
scheme outlined in Fig. 4a. The particle is almost trapped in the resonance sector, 'transverse' 
diffusion due to other resonances is very slow. Time of integration t ~ 6 • 10^. Parameters 
are Ai = —0.02, A2 = 0.05, lo = 4.5 • 10~^. (b) Several low-order resonances are crossed. 
Diffusion in action space is effective due to the zigzag process outlined in Fig. 4b. Parameters 
are: Ai = -0.1, A2 = 0.35, u = 6 ■ 10"^ Time of integration t w 8 • lO'^. 

a way that its trajectory remains almost closed even though geometry of the billiard is 
changing considerably. 

Below we discuss these phenomena in a greater detail. 

Consider the motion of a particle near a (m, n)-resonance. We make a canoni- 
cal transformation (Ji, J2, 0i, ^2) ^ {R,J,(p,4') using the generating function W = 
(m0i — n(f)2)R — {h4'i — h<p2)J 1 where /i, I2 are integers, m/2 — nli = 1. This means 

J = nh + mh, R = hh + hh (6) 

The physical meaning of this transformation is as follows. The new phase = tjk^i — 
n02 (the resonance phase) changes slowly near {m : n)— resonance, while the new phase 
ip = — /i0i + /202 changes fast. New actions R and J are canonically conjugate to and ip, 
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FIG. 6: Diffusion in energy and in action. Left: growth of dispersion of energy of ensemble of 
initially isoenergetic phase points due to crossing of a single resonance as a function of time. 
Right: growth of dispersion of action Ii of the same ensemble. In both cases, values of the 
magnitudes are rescaled to e. Time of integration shown is approximately 100 slow periods, 
with e = 0.0001, a; = 4.5 • 10-^ Ai = -0.02, A2 = 0.05, di = 1, ^2 = 1.3, number of phase 
points is iV = 200. 

respectively. We will see that J is well conserved, while the magnitude of R experiences 
considerable dynamics in the vicinity of the resonance. This fact is useful for the partial 
averaging procedure outlined below. 

The resonance condition now takes the form 

R = Rresi-J, £t) = Tl2~, 2J2~'^ 

171^0.2 + n^df 

In the new variables Hamiltonian is Ti = 'Ho{R,J,et) + uT-Lii^R, J,(j),ip,et). The 
canonically conjugate pairs of variables are {R,(f)) and {J-,ip)- As noted above, it follows 
from the form of the Hamiltonian that ip is a, fast variable, and one can average the 
equations of motion over it. In the averaged system, J is the integral of motion and below 
we consider J as a parameter. We use another canonical transformation (i?, 0) h-> (P, 0) 
with generating function Wi = (P + Rres{J,£t))(f) to introduce new 'action' variable P = 
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R — Rres'i the conjugate 'angle' variable is 0. This transformation allows us to obtain an 
effective Hamiltonian in a vicinity of the resonance, since the new action P measures the 
deviation of R from its resonance value. In a small neighborhood of the resonance where 
P is of order y^, the Hamiltonian takes the following form: 

H = A( J, at) + Fo + 0(£), Fo = ]^gP^ + 6x(0) + (8) 
where A is I-Lq on the resonance, and 



(9) 

8 - 1 / 0-0V7r, if O<0<7r, 

TT^ ^ n"* t + /T'', II — TT < < 0. 

The dynamics near the resonance depends on the properties of the resonance Hamil- 
tonian Fq. Consider this Hamiltonian at frozen values of J,et. Corresponding phase 
portraits of the resonance Hamiltonian on the (0, P)-plane are shown in Figure 1. If 

\a\ < \b\, the phase portrait has a separatrix and an oscillatory domain. If \a\ > \b\, there 
is no such a domain in the phase portrait. To be definite, assume a < 0, 6 > 0. The area 
S of the oscillatory domain (in the case \a\ < \b\) is a function of J and st, and can be 



expressed as follows: 



S = 2J- [ ^FS - bx{(p) - acp d0, (10) 

where Fq is the value of Fq on the separatrix 

TT (g + bf 

^° - 4~^' ^^^^ 

and 01, 02 are the roots of the integrand (in particular, 02 is the coordinate of the saddle 
point). The explicit expression for the area can be found in the Appendix. 
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Now we take into account a slow evolution of parameter et in the system. Suppose the 
area of the oscillation region 5* grows with time. In this case, additional space appears 
inside the oscillatory domain. Hence, phase points can cross the separatrix and change 
their mode of motion from rotational to oscillatory. 

As mentioned above, this phenomenon is called a capture into a resonance. Note 
that for a phase point captured into the resonance, in the course of its motion in the 
phase portrait of Hamiltonian Fq with slowly changing parameters, the internal adiabatic 
invariant (the area encircled by the point's orbit, divided by 27r ) is well conserved (see 
numerical calculations in Ref. {igI). 



It is important to trace the motion of the phase point in the plane of actions (Ji, 12). 

Far from the resonance, a phase point undergoes small oscillations around its adiabatic 
position (a point {l[^\l2^^) in the action plane). 

As time goes, the resonance ray I2 = am;n{£t)Ii slowly sweeps a sector (let us call it the 
resonance sector) on this plane. This resonance ray is a line passing through the origin 
and slowly periodically changing its slope between the minimal and maximal values, thus 
sweeping the resonance sector periodically. When the resonance ray passes through our 
phase point and captures it, the captured point leaves a vicinity of the initial adiabatic 
position R = Rq = hlf"^ + ^1/2°^ = const, J = Jo = const in the action space (Ji, I2) and 
continues its motion following approximately the resonant ray and keeping its value of J 
constant: R = Rres{J,£t), J = Jq = const. Thus, it travels along a resonant trajectory 
which is a part of the line J = Jq located within the resonant sector. 

So, it leaves a vicinity of the initial point /1.2 = -^12 const and continues its motion 
approximately following the corresponding resonant line. Let us denote the value of 
S at the moment of capture as 5*. Somewhen later, the value of 5* starts decreasing 
with time. Approximately at the moment when it equals the value S^, the phase point 
previously captured crosses the separatrix again and leaves the oscillatory domain. This 
is escape from the resonance (value of S at the moment of escape is approximately equal 
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to S** due to conservation of the internal adiabatic invariant). If behaviour of the area 
S{t) is monotonous as the slope am-.nit) is changing between amin and amax, e.g. it 
monotonously grows as am-n is increasing from amin to amax and then monotonously 
decreases decreasing from a^ax and amin, then the escape from the resonance 

happens approximately near the same point (/°, ) where the capture happened. This is 
the case in our numerical examples. 

Phase points that cross the resonance without being captured undergo a jump of adi- 
abatic invariant R of order ^/e (and, therefore, jumps of Ji^2, since Ji = mR — liJ, 
I2 = —nR + I2J). As mentioned above, this phenomenon is a scattering on a resonance. 



Asymptotic formula for the jump in our case is (see also Ref. 2J] for a general formula) 

AR = -2V-e 1 . ^^^.^'^^^'^ , (12) 

where a = a/e, b = b/e, 0* is the phase at which resonance crossing happen, and 
is the value of the Hamiltonian Fq at this crossing (Fq), rescaled to e: = F^/e. The 
value of ^ = Frac[/i*/(27ra^,)] is uniformly distributed on the interval (0, 1) 

Both capture into a resonance and scattering on a resonance are probabilistic phenom- 
ena, and they are usual in systems with resonance crossings (see Ref. |24j. |26|. |27||). 

Multiple crossings through the resonance lead to diffusion of the adiabatic invariants 
due to multiple scattering on the resonance. Note that this diffusion happen along the 
resonance trajectory if only one resonance is taken into account. We discuss this question 
in a more detail in the next Section. 

The phase point undergoes two crossings of the resonance per slow period, and the 
mean value of the change in action is equal, in the first approximation, to the area 5* 
divided by 27r, with a corresponding sign depending on whether the area increases or 
decreases at the moment of crossing. The mean value of the total change in action after 
the two crossings is therefore zero. One can consider the jumps of as a random walk 
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with timestep equal to the half of the slow period and a quasi-random value of the 
is obtained from Eq.f ll2p with assumption of uniform distribution of ^ (Ref. 
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25 



ump 



So, the resulting dynamics in the action space is normal diffusion, at least far from 
the boundaries of the resonant sector sweeped by the resonance ray (near the boundaries, 
motion of the resonance ray slows down and the point effectively "touches" the resonance 
instead of crossing it; numerically, one can see that the point is reflected from the boundary 
of the sector). 

The coefficient of diffusion can be found as follows. If there is no oscillatory domain 
in the phase portrait of Fq, the mean value of the jump f|T2l) is equal to zero, and the 
coefficient of diffusion is equal to the expectation value of (Ai?)^, divided by half of the 
slow period, similar to that being done in Ref. [26[. It can be shown that the quasirandom 
variable ^ is uniformly distributed on (0, 1) and its values on consequent passages through 
the resonance are independent. In case there is the oscillatory domain in the phase portrait 
of Fq, the mean value of the jump AR is not equal to zero, and the coefficient of diffusion 
can be found by calculating the expectation value of (AR^^^ + AR^'^^y from (fT2|) (where 
AR^^\ AR^'^^ are jumps of R at two consequent passages through the resonance, with 
values of ^ equal to .^i and ^2, correspondingly, assuming the same uniform distribution 
of the quasirandom variables ^1 and ^2), and dividing it by the slow period. 

Numerical calculations are presented in Figs. 2-3. For numerical simulations we chose 
parameters di harmonically varying with time. In Figure 2 jumps of adiabatic invariant Ji 
are shown. A single jump of the adiabatic invariant is shown in the inset of Figure 2. In 
Figure 3a one can see the evolution of the adiabatic invariant due to multiple scatterings 
on the resonance and one capture of the phase point into the resonance. The captured 
point moves along a resonant curve until it escapes from the resonance. Figure 3b gives 
evolution of Ji, I2 in the vicinity of capture into resonance. Fig. 3c presents time evolution 
of the winding number bJxjbJi- 
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IV. ACCELERATION BY SEVERAL RESONANCES 

Consider what is happening if parameters di{t) are changed periodically in such a way 
that the particle repeatedly cross a single resonance. I.e., only one low-order resonance, 
say n:m, is crossed, while other resonances being crossed are of high order and can be 
neglected (While there is no precise distinction between low- and high-order resonances, 
in this paper we consider as high order resonances those for which |m| + |n| is bigger than 
const-e"^''^. This is a condition that the resonance alone does not produce any diffusion 
because its effect is so small that phases of subsequent passages through this resonance 



are strongly correlated, see also Ref. [26|). Periodic scattering on the low-order resonance 
leads to diffusion in action space. However, jumps of Ji and I2 at the passage through 
the resonance are linearly dependent. Indeed, since J = nli + ml2 ~const at the passage 
through the resonance, jumps of Ji and I2 are dependent in such a way that in the action 
plane (hyh) the particle is kicked along the corresponding line, J = Jq =const. Since 
the resonance condition is I2 = am:n{£t)Ii, one can easily see that the resonance ray will 
periodically sweep a sector on the plane of actions, and the particle will diffuse only along 
the part of the line J = nli + m/2 =const restricted by this sector. In other words, 
its diffusion in action space happens along the interval. This is an obvious geometric 
obstacle to unlimited acceleration. Of course, taking into account that near a low-order 
resonance there are always high-order resonances, one gets diffusion transverse to the 
interval J = Jq, although very slow. Now, to let the energy of the particle to grow 
without bounds efficiently, we should allow at least two low-order resonances (say, n, m 
and n, m), and let the corresponding resonance sectors to overlap. Then, in the first sector 
diffusion happens along the line nli + ml2 =const and in the second sector along the line 
nil + mJ2 =const. Zigzag- type acceleration is possible: shift of the particle in the action 
space due to the first resonance is not parallel to the shift of the particle due to the other 
resonance. Even though the particle is trapped within the union of the two resonance 
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sectors, it is now can drift to infinity staying within this combined sector. Schematically 
this mechanism is presented in Fig.4. 

Numerical examples are given in Figs 5-6. In Figure 5a one can see diffusion of the 
phase point in the action space restricted by a resonance sector. The particle is effectively 
trapped on an interval, and diffusion transverse to this interval is extremely slow. In Fig. 
5b, due to subsequent crossings of several resonances, the particle can undergo unlimited 
diffusion. From stroboscopic shots of the position of the particle in action space shown 
in this Figure, peculiar properties of the motion of the particle can be recognized. There 
are regimes of motion where the particle diffuses along certain direction for a long time 
(when it is in a region of the action space swept by only one low order resonance), and 
there are regimes of motion, where it receives subsequent kicks in different directions (in 
the region of overlapping of resonance sectors). 

Fig. 6a demonstrates growth of dispersion in energy of an ensemble of phase points 
due to multiple scattering on a single resonance. It is seen that the initial growth is linear, 
in accordance with the theory for the single-resonance system. However, saturation or 
only extremely slow diffusion is expected after the ensemble of particles is redistributed 
within the restricted interval shown in Fig. 5a. 

Similar results can be seen in Fig. 6b, where multiple scattering on three low-order 
resonances happens. Since diffusion of energy is bounded from below, we expect the 
total energy of the ensemble will also start grow with time after considerable spreading 
in energy happens. 

It should be noted that the reasoning in this Section is of heuristic nature. While for the 
partially averaged system the formula Eq.f ll2l) for the jump of action at the resonance and 
the statement about conservation of J are mathematically rigorous results, the using of 
these results in the exact two-frequency system is based on physically plausible reasoning. 
We use numerical simulation in order to support heuristic reasoning. 
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V. CONCLUSION 



An important open question left in Ref. [16| concerned the possibihty of unhmited 
acceleration of a particle in the billiard under consideration. It is well known (see, e.g. 
Ref. l"^) that in a similar one- dimensional problem (Ulam's model (Ref. motion of 
a particle between two periodically moving walls) unlimited acceleration is impossible, 
provided that the motion of the walls is described by smooth enough functions. The 
reason is that when a particle moves fast enough with respect to the motion of the walls, 
the system possesses a perpetual adiabatic invariant (see Ref. {2^) that prevents the 
particle's acceleration. In the two-dimensional model considered in the present paper 
and in Ref. 16|, the adiabatic invariance is destroyed by captures into resonances and 
scattering on resonances. The hypothesis of Ref. [l^ is that for a majority of initial 
conditions the velocity of a particle can reach arbitrarily large magnitudes. However, 
there was not clear numerical evidence for that in Ref. [16|]. 

Here we show that there are geometric obstacles to unlimited acceleration in the case 
where only a single resonance influences dynamics of the particle. The particle is then 
trapped in the resonance sector. Overcoming this obstacle by modulating the system in 
such a way that several low-order resonances are crossed during a slow period of modula- 
tion, one obtains clear numerical evidence of unlimited acceleration. 
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Appendix 



A. Fourier series 



Variables pi and from Eql2] can be expressed in terms of (/j, (pi) as the following 
Fourier series: 



P^ = TTsinA;0i, (13) 

k=l ^i'^ 



.=1 ^'k^ 



where primes denote summation over odd k. 

Note that U is related to pi in a very simple way: 

, '^d\pi\ 



(14) 



as can easily be seen from the phase portrait of the one-dimensional billiard. From this 
expression it is easy to obtain y = which give us the first term of Eq.(I3]). 



B. Area of the oscillatory domain 

The explicit expression for (pi in Eq. lITUi) depends on the relation between a and h. 
Consider, e.g. the case a < {2\/2 — 3)6. One has 

Tx a + h r- Ti a + h 

01 = -2-^(1 + V2)> '^2=2-^- (15) 

For the area 5" one finds: 

S=,/I(^fl + ^Un'^ (16) 



g 2 V ij\bj 
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